Survival of a Diffusing Particle in a Transverse Shear Flow: A First-Passage Problem 
with Continuously Varying Persistence Exponent 
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We consider a particle diffusing in the indirection, dy/dt — n(t), subject to a transverse shear 
flow in the ^-direction, dx/dt = f(y), where x > and x = is an absorbing boundary. We treat the 
class of models defined by f(y) = ±v±(±y) a where the upper (lower) sign refers to y > (y < 0). 
We show that the particle survives with probability Q(t) ~ t~ e with 9 = 1/4, independent of a, if 
n+ = V—. If v+ 7^ however, we show that 8 depends on both a and the ratio v+/v~, and we 
determine this dependence. 

PACS numbers: 02.50.-r, 05.40.-a 
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First-passage problems are an important aspect of the 
theory of stochastic processes. Applications include pop- 
ulation dynamics, chemical reactions, and many other 
problems in science [1]. This class of problems has at- 
tracted a resurgence of interest in the last decade. The 
number of exactly-solved models, however, is still quite 
small [1,2]. The classic example is a one-dimensional 
random walker, dx/dt = f](t), where rj(t) is a Gaussian 
white noise with mean zero and correlator (f](t)r](t')) = 
2DS(t — t'). If the particle starts at x > 0, and there is 
an absorbing boundary at x = 0, the particle survives at 
time t with probability Q{x,t) = erf (a; / \/4Dt) [1]. For 
large t this gives Q(x,t) ~ x/(Dty/ 2 , defining an expo- 
nent 9 = 1/2 (often called a 'persistence exponent' [2]) 
through the relation Q(x, t) ~ t~ e for large t. 

In this paper we introduce a new class of models for 
which 9 can be determined exactly and has, in general, 
a nontrivial value. We treat a simple model of a parti- 
cle diffusing in the ^-direction, dy/dt = 7]{t), but sub- 
ject to a transverse flow field ('shear') in the cc-direction, 
dx/dt = f(y), with an absorbing boundary at x = 0. We 
show that, in general, the exponent 9 depends continu- 
ously on the model parameters, except when the function 
f(y) is an odd function, f(—y) = —f(y), when we find 
9 = 1/4. Two specific examples with odd f(y) have 
attracted some attention in the literature. The first, 
f(y) = Z/j is equivalent to d 2 x/dt 2 = rj(t), the 'random 
acceleration model', for which the survival probability is 
known rigorously to decay as Q(t) ~ t -1 / 4 [3]. The sec- 
ond model for which we are aware of previous work is the 
case f(y) = v sgn(y). Redner and Krapivsky [4] have pre- 
sented qualitative arguments that 9=1/4 for this model 
too, and indeed have conjectured that 9 = 1/4 for all odd 
functions f{y). Below we present an argument, based on 
the Sparre Andersen theorem [5], that this conjecture is 
correct. First, however, we derive the result analytically 
for the class of odd functions f(y) = vsgn(y)\y\ a . Fur- 
thermore, we will show that when the amplitude v takes 
different values, u+ and u_, for y > and y < re- 
spectively, the exponent 9 becomes a nontrivial function, 
whose form we determine explicitly, of the exponent a 



and the ratio w+/u_ . Our result thus greatly extends 
the class of first-passage problems for which the value of 
9 is known exactly. 

The starting point of the calculation is the two equa- 
tions 



y = v(t) , 

X = f(y) = v±sgn(y) \y\ a 



(1) 

(2) 



where the upper (lower) sign refers to y > (y < 0), 
dots indicate time derivatives, and rj(t) is Gaussian white 
noise. The final result for 9 can be written in the form 
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where j3 = 1/(2 + a). Note that when v + = w_ the result 
reduces to 9 = 1/4, independent of a, as claimed. The 
general result is easily understood in the limits i>_ — > 
and v+ — > 0. For t>_ — > 0, there is no mechanism by 
which the particle can reach the absorbing boundary at 
x = 0, so 9 = as predicted by (3). For v + — ► 0, in order 
to survive the particle should avoid the y < region. The 
probability to remain in the y > region decays as i -1 / 2 
as in the one-dimensional random walk, i.e. 9 = 1/2, 
again in accord with (3). 

We now outline the derivation of Eq. (3). FromEqs. (1) 
and (2) we can write down the backward Fokker-Planck 
equation 



dQ 
dt 



= D 



d 2 Q 
dy 2 



± v±(±y) 



dx 



(4) 



where Q(x,y,t) is the probability that the particle still 
survives at time t given that it started at position (x, 
y), and we will take x > 0. The absorbing boundary at 
x = leads to the boundary condition Q(Q,y,t) = for 
all y < 0. The initial condition is Q(x,y,0) = 1 for all 
x > 0. 

Solving the full initial value problem is difficult so 
we will specialize to the late-time scaling regime where 
Q(x,y,t) ~ t~ e . This approach exploits a generaliza- 
tion of the method recently introduced by Burkhardt for 
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a = 1 (the random acceleration problem) [6]. Dimen- 
sional arguments give us the appropriate variables for the 
problem. From (1) and (2) we obtain y ~ (Dt) 1 / 2 and 



v(Dt) 



a/2 



giving x 



(v/D)y 1 ^ 3 , where (3 = l/(2 + a) 



and v carries the dimension of the amplitudes v±. The 
combination z = vy 1 /^ /Dx serves as a dimensionless 
scaling variable, which can also be seen directly from Eq. 
(4). This gives, for large t, Q(x,t,y) ~ [y 2 /Dt) e G{z) 
where G(z) is a scaling function. This is a little loose, 
however, as we expect different solutions for y > and 
y < due to the non-analyticity [7] of f(y) at y = 0. 
Taking this into account, and expressing the prefactor 
of the scaling function in terms of x rather than y for 
convenience, we write, for t — > oo, 



Q(x,y,t) 



,.23 



F±[± 



v±(±y) 



V/3 



Dx 



(5) 



where F±(z) is the scaling function for y > (+) or 
y < (— ). The (dimensional) prefactors (for y > and 
y < 0) in Eq. (5) have been omitted since Eq. (4) is lin- 
ear. The functions F+(z) and F-(z) are defined such 
that the prefactor is the same for y > and y < 0. 

Inserting the form (5) into the backward Fokker-Planck 
equation (4), we see immediately that the term dQ/dt 
leads to a term of order t~( B+1 \ which is subdominant 
for large t and can therefore be dropped. The remaining 
terms give 

zF'l{z) + (1-/3- l3 2 z)F' ± (z) + 2f3 3 8F±(z) = . (6) 

Expressed in terms of the variable u — (3 2 z, this equa- 
tion becomes Kummer's equation. Independent solutions 
are the confluent hypergeometric functions M(— 2/30, 1 — 
P, /3 2 z) and 17 (-2/30, 1 - /3, /3 2 z) [8]. The appropriate so- 
lutions for each domain are selected from the limiting 
behavior as x — > 0. For y > 0, we require Q(0,y,t) <~ 
(y 2 /t) d which implies, from Eq. (5), that F+(z) ~ z 2f39 
for large z. The function M(— 2/30, 1 — (3,(3 2 z) diverges 
exponentially for z — ► oo, so must be rejected for 
y > 0. For y < 0, Q must vanish at x = 0+ (be- 
cause the flow field immediately takes the particle on 
to the absorbing boundary). For large negative z, both 
M(-2/30, 1 - /3, (3 2 z) and U{-2/39, 1 - /3, /3 2 z) behave as 
(— z) 2/3e , so in this regime we retain both solutions and 
fix 9 by the requirement that the coefficient of (— z) 2/3e 
vanishes. Thus we write 



F+(z) = AU \-206, 1-/3, 
F_(*) = BU (-2/36,1- /3, 



+CM [-2/30, 1-/3,- 



Dx 
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(8) 



Relations between the coefficients A, 7?, and C can 
be obtained by imposing the appropriate continuity con- 
ditions at y = 0. Requiring F+(Q) = 7-(0) gives, 



using M(o, b, 0) = 1 and, for b < 1, ?7(o,6,0) = 
7r[sin(7r6)r(6)r(l + a- fe)]- 1 [8], 



C = (A - B) 



sin(7r/3) r[l - /3]r[/3(l - 26*)] ' 



(9) 



for all /3 > 0. The second relation between A, B, and C 
is obtained by noting that, for a > — 1, Eq. (4) implies 
that the first derivative, dQ/dy, is continuous at y = 0. 
The term linear in y comes from a term of order z 1 ~ b in 
the the small- z expansion of U (a, b, z). Since 6=1-/3 
here, this term is of order z 13 , i.e. linear in y since the 
scaling variable z is proportional to y 1 ^ , as in Eq. (5). 
Demanding that this linear term has the same coefficient 
in both regimes (y > and y < 0) yields 



Av'l = Bit 



(10) 



Finally we need the asymptotic behavior of F_ (z) for 
z — > — oo. From the asymptotic properties of the func- 
tions M(a, b, z) and ?7(o, 6, z) [8] one readily obtains 



F-{z) 



B 

x(- 



[sin[7r/3(l-40)/2] 
sin[7r/3/2] 



r[i - /3] 



r[i-/3(i 



— oo. 
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Since F_(z) should vanish for z — > — oo, due to the ab- 
sorbing boundary at x = 0, the coefficient of (— 2;) 2/3e 
must vanish in Eq. (11). This condition, combined with 
Eqs. (9) and (10), gives the result presented in Eq. (3). 

The case v + = where f(y) is an odd function, is 
especially simple. Eqs. (10) and (9) give A = B and 
C = respectively for this case whence, from Eq. (11), 
one obtains sm[nf3(l — 40)/2] = 0. The desired solution 
corresponds to the smallest positive value for 9, which is 
= 1/4. Redncr and Krapivsky [4] have conjectured that 
9=1/4 for all odd functions f(y). A general argument 
for this result can be made using the Sparre Andersen 
theorem [5]. The theorem deals with a one-dimensional 
random walk, with an absorbing boundary at x = and 
step sizes that are independent, identically distributed 
random variables drawn from any continuous, symmetric 
distribution. The result of the theorem is that the proba- 
bility that the walker still survives after N steps decreases 
as iV -1 / 2 for large N. To apply this result to our prob- 
lem, we treat the successive crossings of the x-axis by our 
diffusing particle as defining the steps of the walk, i.e. 
successive crossings at x = x n and x = x n+ \ correspond 
to a step of length x n +i — x n in an effective random walk 
along the x-axis. Clearly the distribution of step sizes 
is continuous, and is symmetric for any odd f(y). Since 
the number of crossings in time t scales as N ~ t 1 / 2 , the 
survival probability decays as Q ~ N~ x / 2 ~ i -1 / 4 [9]. If 
v + ^ V- , the distribution of step sizes is no longer sym- 
metric, the conditions of the Sparre Andersen theorem 
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do not hold, and 9 is different from 1/4. A superficially 
related model is the Matheron-de Marsily (MdM) model 
[10], introduced to model the dispersion of tracer parti- 
cles in porous rocks. In this model, the function /(y) in 
Eq. (2) is a random function of y with zero mean. The 
MdM model also has 6 = 1/4 [11], but the Sparre An- 
derson theorem cannot naively be applied because /(y) 
is not odd. 

Returning to our original model, we note that two spe- 
cial values of a are of particular interest. The case a = 
corresponds to the original 'windy cliff' model of Rcdner 
and Krapivsky [4], generalized to different drift velocities, 
v + and V-, for y positive and y negative. If we transform 
to a frame of reference moving at speed u along the x- 
axis, then in the new frame the flow velocities have values 
±v± — u, while the absorbing boundary has velocity — u. 
The expression for 9 presented in Eq. (3) therefore applies 
equally well to this class of moving boundary problems. 
The second interesting value is the case a — 1. Elim- 
inating y from Eqs. (1) and (2), one obtains a random 
acceleration process x = v±r)(t), with different strength 
noise forces depending on the sign of the velocity. 
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FIG. 1. Simulation data for 9 for a range of values of the 
exponent a and amplitude ratio v+/v- in Eq. (2). The con- 
tinuous curves are the predictions of Eq. (3). The errors on 
the data points are smaller than the sizes of the symbols. 

The result (3) has been checked by computer simula- 
tions for a range of a and The results are pre- 
sented in Fig.l. Motion in the y-direction was modeled 
by a discrete-time lattice random walk, while the x co- 
ordinate is treated as a continuous variable. The initial 
condition was x = 0, y = 0, with the first step taken in 
the positive y direction. The simulations enter the power- 
law regime, where Q(t) ~ t~ e , quite quickly, making it 
relatively straightforward to obtain an accurate estimate 
of the exponent 9. The data (symbols) are in excellent 



agreement with the predictions (lines) in all cases. Notice 
that all curves and data sets cross at 9 = 1/4, v+/v- = 1, 
as anticipated. 

The present exact approach is restricted to the case 
where the flow velocity f(y) has a pure power-law form, 
with the same power, for both y < and y > 0, although 
the amplitudes in the the two regimes can be different. 
What can we say about the general case? The case of 
different powers for y > and y < cannot be treated 
analytically within the present approach because the as- 
sumed scaling form, Eq. (5), for Q(x, y, t) does not lead to 
a consistent solution for both signs of y. Consider, how- 
ever, the case where the function f(y) has, for y — ► ±oo, 
an asymptotic power-law form with the same power for 
both signs of y. Since trajectories that survive for a long 
time explore large values of \y\ [4] we conjecture that the 
exponent 9 describing the asymptotic decay of Q{t) is 
given by Eq. (3) with the exponent a and amplitude ra- 
tio tL|_/i>_ extracted from the limiting behavior of /(y) 
for y — ► ±oo. This line of reasoning leads immediately to 
the further conjecture that when the asymptotic behav- 
ior of /(y) for y — ► ±oo is described by different powers, 
a±, the larger power dominates the t —* oo behavior, 
such that if a + > a_ then 9 — (i.e. there is a non-zero 
probability that the particle escapes to infinity), while if 
a + < a- then 9 = 1/2. 

One may also consider the effect of diffusion in the in- 
direction as well as the y-direction. For the case f(y) = 
v sgn(y), corresponding to a = 0, it was argued in ref. [4] 
that the diffusion in the x-direction is asymptotically ir- 
relevant compared to drift since Xdis ~ (Dt) 1 ^ 2 whereas 
Adrift ~ vt. Applying this reasoning to general a gives 
Adrift ~ i( 2+Q )/ 2 , so drift in the ^-direction dominates 
diffusion for all a > — 1, the latter condition coinciding 
with the range of a for which our analytic approach is 
valid. We conclude that diffusion in the ^-direction does 
not change the value of 9 (for a > — 1). A more formal 
way to see this is to introduce a term D x d 2 Q/dx 2 , corre- 
sponding to diffusion in the x-direction, on the right-hand 
side of Eq. (4). Under the scale change y — > Ay, t — > X 2 t, 
and x — > X 1 '^x, the equation retains its form except for 
a rescaled value of D x , given by D' x = A~ 2 ( 1+Q ) D x . This 
shows that D x scales to zero on large length and time 
scales and is therefore asymptotically irrelevant. 

In summary, we have obtained exact results for the 
persistence exponent 9 for a class of two-dimensional 
stochastic processes in which a particle diffuses in the 
y-dircction and has a deterministic, y-dependent drift in 
the x-direction and an absorbing boundary at a; = 0. 
The universal exponent 9 = 1/4 is obtained when the 
drift velocity is an odd function of y, but 9 has a non- 
trivial value, between zero and 1/2, in other cases. These 
results significantly expand the (up till now) rather small 
class of soluble first-passage problems with non-trivial ex- 
ponents. 

AJB thanks Richard Blythe for useful discussions. 
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